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ABSTRACT 



A need for more knowledge of meteorological parameters in 
the region 30,000 to 45,000 feet has cpme about in recent years 
as a result of rapidly growing jet aircraft operations. This 
paper discusses one of these parameters, the height of the 
level of maximum winds, which affords considerable aid to the 
flight planner in his selection of preferred flight altitudes 
and routes. The deviations of the level of maximum winds from 
the tropopause are discussed and graphically shown, indicating 
that the latter is not a satisfactory first approximation to 
height of the former. This fact gives rise to the development 
of a series of simple and multiple regression equations for 
determining the height of the level of maximum winds as a 
function of both the tropopause and the 300-millibar heights. 

Professor Robert J. Renard of the U. S. Naval Post- 
graduate School has been the guiding force behind this investi- 
gation. For his time and effort, the writer is deeply 
indebted. 
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1. Introduction. 



With the advent of jet aircraft in recent years, the fore- 
caster’s spatial envelope of responsibility has almost doubled. 

It has become necessary to analyze the meteorological parameters 
up to 100 mb on a daily basis. Knowledge of v/ind, cloud, tem- 
perature, and pressure fields at these levels are necessitated 
by the demands of aircraft flight planning. The United States 
Weather Bureau [^9^ has suggested that the forecast accuracy of 
the temperature and wind at jet flight levels be + and 
+ 20 knots respectively. 

Not the least important of the parameters needed for econom- 
ical operation of aircraft over long distances are the v/ind 
speed and direction, including their variation both in time and 
space. Over the continental United States and parts of Canada, 
a United States Weather Bureau analysis of the layer of maximum 
wind (hereafter abbreviated LMW and used synonymously u'ith level 
of maximum wind) is available on a 1/40,000,000 polar stereo- 
graphic facsimile chart every 12 hours. The height of the UfW 
typically varies from about 24,000 feet in the northern sec- 
tions “to about 44,000 feet in the south. This analysis, along 
with the contour analyses of the standard isobaric surfaces, 
reasonably describes the horizontal wind field and its deriva- 
tives throughout the high troposphere and lower synoptic strato- 
sphere,^ However, the extension of the LMW analysis, particularly 

a term defined by H. A, Panofsky [lOj , which includes the 

region from the tropopause to about 10 millibars. 
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to ocean areas, can be accomplished only with great difficulty 
owing to the paucity of wind data in the stratosphere. Thus, 
some means of determining with reasonable accuracy, the LMW data 
from the more easily obtained variables at lower levels would 
allow extension of the high-level wind analysis to ocean areas. 
One phase of the problem involves determination of the height 
of the LMl^. It is this parameter which is under statistical 
investigation here. 
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2. Background, 

Much has been written concerning the eiivironiTient of tiie jet 
aircraft since the inauguration of high altitude, long range 
flights within the past ten years. Relative to the specific 
problem at hand, Reiter ^5^ discusses many aspects of the layer 
of maximum winds, including reliability of upper-wind data, 
analysis of the horizontal and vertical structure of the LMV7, 
and forecasting the LMW, with special application to jet flight 
operations. His general conclusions testify to the need for not 
only more accuracy in upper -wind measuring devices, but to their 
proper representation and interpretation as well. 

Since there are rather large inaccuracies in the upper- 
wind data and because of the scarcity of data in certain areas, 
it becomes desirable to determine the LMW from more commonly avail- 
able data. The United States Weather Bureau assumes that, for 
purposes of analysis, the height of the LMW is at, or very nearly 
at, the tropopause in middle and northern latitudes. The tropo- 
pause, if considered a first-order discontinuity in the temper«a- 
ture field, is logically where one v;ould expect to find the 
maximum \>rind. Unfortunately, as will be shown later, there are 
many occasions in which the LMW deviates significantly from the 
tropopause height, Johannessen C^J discusses a series of ob- 
servations taken at Larkhill, England, v/hich indicates that 
the LMT'J occurs from 20Q0 to 4000 feet below the tropopause. 

However, extension to other latitudes is not plausible. 

The investigative results whose discussion follows is a 
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consequence of the actual deviations of the height of the 
from that of the tropopause over the United States during the 
months of November 1959 and January I960, which include a 
great variety of weather types. 
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3, Deviations of the Height of the LMW from the Height of the 
Tropopause 

The tropopause definition as used by the United States Weather 
Bureau [iQ is, using only the lapse-rate criteria, that point in 
the sounding at which the lapse rate decreases to 2°C/km and then 
averages less than this in the 2 -km layer immediately above. 

Whether nationally or internationally defined, the tropopause 
level is a controversial subject. Stinson {s] in recently proposing 
a new tropopause identification and classification system deals 
with a myriad of soundings which bring out the difficulties in 
selecting a certain point fulfilling the above lapse-rate condition. 

Regardless of inaccuracies or uncertainties of tropopause 
location, the heights as given on the United States Weather Bureau 
OOZ and 12Z facsimile charts of tropopause heights (Z^) were used. 
These charts were analyzed by the author for use in this investiga- 
tion. A typical example of the analysis is shown as part of 
Figure 1. A striking feature is the band of strong height gradient 
oriented generally east-west from southern California through 
Texas to Florida. This represents, for the purpose of this study, 
a continuous zone of transition between the tropical and middle - 
latitude tropopauses. Undulations of this band are always present 
in varying degrees of amplitude and are associated with the major 
weather systems. For example, the passage of a low or trough at 
300 mb is synonymous with a trough in the tropopause heights. 

Another zone of large height gradient appears in northwest United 
States and is probably associated with the transition from the 
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iaiddic to the polar tropupause* A hemispheric analysis oi 
this same phenomenon [3^ shown in Figure 2, clearly exhibits 
Che nature of Che two transition scones just mentioned* 

As an approach Co defining statistically some relations 
betv/een the height of tlie 1*MV/ and the tropopausu, the 

United States Weather Bureau analysis of the height of the LMW 
(see Figure 1) is compared to the locally produced tropopausc 
analysis by a graphical subtraction, of the contour heights; of 
tlie two surfaces. This also is shown in Figure 1. The general 
area of approximate coincidence of the tv;o surfaces extends 
from- extreme northwest United States tov/ard Che southeast 
through Missouri and from extrevae northeast United States 
toward the southwest through Tennessee. In general, this area 
of coincidence, defined as that section in which 
follows closely Che 36,000 to 40-,000-ft channel in the Cropo- 
pause field* In areas of** tropopause height less than 36,000 
feet, while tropopause heights* greater than 

40,000 feet are associated with Further, since 

the tropopause is closely associated with major weather systems, 
it is reasonable to suspect that the* troughs, ridges, and 
centers in Figure 1 will have continuity from day to day. That 
this is true is shown in Figure 3,’ which is a graphical sub- 
traction for the situation as it existed 36 hours previous 
to that in Figure 1, The movement of the minus eight isoline 
over Montana corresponds closely to the ’movement of a 
27,600-ft lov; at Che 300-mb surface situated about 300 miles 
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Fig. 2 Tropopause Height Contours Jan. 1 , 1956 . 

(After Defant and Taba.) 
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to the north. The plus 16 area moved east and increased to a 
plus 20 over Texas with another plus 16 coming in over central 
California on Nov 17. These large plus areas appear to follow 
the movement of a 300-mb ridge located 400 miles to the west. 
Thus, figures 1 and 3 show that an assumption of equalling 
made only over a limited area of the continental 
United States. For ^the two months considered, the tropopause 
was found to vary from 12,000 feet below the LMW in northern 
United States to as high as 24,000 feet above the LMW in the 
south. ^ ‘ . 
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4, Statistical Relationship Between the Height of the L^^/J and 
other Parameters 

A study of Figure 1 suggests a definite relationship between 
(Zt“Zlj^j) and Zj. A better picture of this relationship is 
shovm in Figure 4 as a scatter diagram of Zj versus 
the points being randomly selected from November 1959 and January 
1960. A best-fit line may be drawn through these points using 
the least squares method but because of the decrease in slope for 
the envelope of data below = 36,000 feet, the standard error 
would undoubtedly be excessive. Thus, two straight lines may be 
drawn, one for Zj greater than 36,000 feet and the other for Z^^ 
less than 36,000 feet. The discontinuity which occurs betv/een 
the two regression lines at Z^ = 36,000 feet can easily be smoothed 
in to provide one continuous curve. Before proceeding on computa- 
tional v;ork, approximately 1000 actual data points, distributed 
quite evenly over the two-month period, were selected, and data 
for Zj versus Zj 2 j\^ recorded. About half of these randomly- 
selected points v;ere plotted on the scatter diagram of Figure 5. 
This diagram shows that a two-dimensional linear relationship 
will give best results for values of Zx ranging from 30,000 feet 
to 44,000 feet. Above and below this range the slope of the' 
data envelope changes and the scatter increases. 

j 

With a view toward simple utilization of all data points, 
it was felt that the best course of action was to work first with 
Che data as originally divided, that is, for tropopausc heights 
above and below 36,000 feet. In order to use the 1000 points 
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taken over the period ^ it became necessary to employ the 
National Cash Register 102A electronic digital computer at the 
United States Naval Postgraduate School. The computer was 
programmed to accept a maximum sample size of 256. Outputs of 
the mean, standard deviation, variance, and linear correlations 
between the tV7o sets of data were recorded. The correlation 
coefficient (r) used is defined in Hoel [^ij as 



p — < ~x ) iYi -'P) 

where s^ and Sy are* the standard deviations of the two varia- 
bles involved. The first correlation attempted related the 
independent variable, Z-j less than 36,000 feet, to the dependent 
variable, (Zj - henceforth referred to as 4 Z. A 

sample size of 256 randomly selected points was used.^ This 
group of data was part of the 1000 points previously mentioned. 

A correlation coefficient of 0.66 resulted. The form of the 
regression equation, utilizing r is 

- -h h 
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The resulting equation is 



f8~S~ St) 

Solving for Zlmw 

The standard error is 2490 feet and can be interpreted as the 
^LMW distance above and below the regression line which includes 
about 687. of the data points. The standard error defined by 
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Panotsky is siven by the expression 

(^i I 

During the course of the graphical subtraction mentioned 
earlier^ it v;as noted that there exists a marked resemblance 
between the contours of the LMW and the 300-mb chart. Thus, for 
each of the 1000 points over the two months* period, a 300-mb 
height was determined and plotted against on the scatter 

diagram shown in Figure 6. The relation here is about the 
same as that for Zj versus ^ Z. However, the 300-mb versus 
Zx plot does not show the Vadical change in slope associated 
with the scatter diagram of Zj plotted against Using 

only those points associated with less than 36,000 feet, 
the correlation coefficient, ^lMW, 300, 0.67. The regression 

equation, computed as before, is 

^ ^ j y/7 /ePO's (Z) 

with a standard error of 2230 feet. The graphs of equations 

(1) and (2) are shown in Figure 7, 

If the t^^’'o independent variables Z-p and Z^qq do not account 
for the same fraction of the variance of Zy^^, there will he 
some improvement in accuracy by writing a multiple regression 
equation of the form 

(PL i-' c ^'3) 

where a, b, and c arc constants to be determined from data in 

the sample. This best-fit plane may be determined by the 

least squares method discussed in Panofsky hy subtracting 
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Figure 7. Graph of simple linear equations applicable to 



the observed value of the height of the LMW, from each 



LMW^ 



side of the above equation, and squaring as 




Partially differentiating the equation with respect to a, b, and 
c successively, gives three equations v/ith a, b, and c as the 
three unknowns o All other parameters in the equation are knov/n 
or can be calculated from the previously mentioned results 
given by the electronic computer^ The simultaneous solution of 
these three equations in this particular case is a very arduous 
task requiring that an excessive number of digits be carried 
along in order to arrive at a reasonably accurate solution. 

Since the number of digits required exceeds the keyboard capacity 
of the ordinary Marchant desk calculator, the symbolic solution 
for constants b and jC as given by Panofsky [aJ is more practical 
and equally as valid. These equations are as follows: 



The constant a may Chen be found from the following expression: 



Thus, for a solution, the simple linear correlation coefficients 
betv/een all the variables taken two at a time, the standard 
deviation, and the means of all the variables, must be known. 
This method makes maximum use of the work previously performed 
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by the coraputei* and eliminates many hours of error-prone v/ork 
on the desk calculator. The computer results used in calculating 
this multiple regression equation are shown in Tables 1 and 2. 





A Z 




^300 


Zt 


0 660 






^300 


.061 


.565 






^LMW 


— 


.580 


.671 



Table 1, Correlation Coefficients for 





Mean 


S tandard 
Deviation 


AZ 


-1,215 


3,313 


^LMW 


32,338 


3,009 


^300 


29,363 


758 




31,064 


3,932 



Table 2, Means and Standard Deviations of 
all variables ior iSj- JOyOaa' 



Thus the regression equation (3), applicable for Z,j, less 
than 36,000 feet, becomes 

'7 ,rf /ao^ ft) 

Here the multiple correlation coefficient, R 
is 0.713 and the standard error is 2110 feet. 

Using 4Z as the dependent variable and solving the 

resulting equation for gives 

2r- i^^/Ciaa^es i^t) (t) 
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v/hich has a 0^760 but shov/s no improvement 

in the standard error because of the larger variance of the 

I' ■ 

dependent variable Nomograms for equations (6) and (7) 

are shown in Figures 8 and 9. 

The multiple correlation coefficient via Panofsky &.] 
is the square root of the following equation: 

. 3oo - ^ 

i- 

The quantity gives that fraction of the total variance of 

the dependent variable that is accounted for by the two 

■» 

independent variables. 

The standard error for a three-variable linear regression 
equation was calculated by a method given by Snedecor 
and involves the use of the following equation: 






: P ^ r LLr.P.n A0 
" L 









v;here n is the sample size and n-3 represents the number of 
degrees of freedom. If the sample size is large, the quantity 
utAsP fj.ri»v)niay be accurately represented by 

^-3 n 

the variance of the dependent variable. The equation reduces 



simply to 



s 



cJ 






The quantity(l - r 2 ) gives that fraction of the variance that 
is not accounted for by the two independent variables and is 
commonly referred to as the coefficient of alienation. 

The section of the data for which is greater than 
36,000 feet has been treated in exactly the same manner as 
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Figure 9. Nomogram for = —5'J2 3 'S'oS Er r ->--3--37 ^ 



herecolore described for Zj less than 36,000 feet. Tables 3 
and 4 summarize the electronic computer results for a sample 
size of 256. 





AZ 


Zrj. 

1 '■ 


^300 




.845 






^300 


.196 


.555 




^-LMW 




.336 


• 662 



Table 3. Correlation Coefficients for ooa^ 





Mean 


Standard 

Deviation 


d z 


4,295 


5,611 


^LMW 


38,504 


3,199 


^300 


30,611 


600 


Z™ 


42,771 


5,865 



Table 4„ Means and Standard Deviations of 
all Variables for SOjOoo'’ 



The equations ca culated from the data of Tables 3 and 4 
and their standard errors are as follows: 



~ '5-t . - -3CJ/0 i=r 

vS-/:r = 

•= ~ 5,£L -- 

rv 9/6> ■h'S'Bt 



i8) 

(9) 

(/o) 

(//) 
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Equations (8) and (11) were obtained using as the dependent 
variable o The graphs of equations (8) and (9) have been dra^^i 
in Figure 10 » Nomograms have been prepared for the multiple 
regression equations and are displayed in Figures 11 and 12, 

One additional effort was made to obtain regression equa- 
tions with smaller standard errors. This was attempted, using 
the same independent variables, by choosing data over ranges 
of values of the independent variables which have the most 
apparent linear relationship with the dependent variable. The 
points randomly selected foV this correlation were taken from 
the original 1000 points and were those for which 

• Figures 5 and 6 

shov7 these ranges of values to be a reasonable selection. The 
computer results using one sample of size 256 are summarized 
in Tables 5 and 6 below: 







^300 


^'LMW 


.694 


.696 


^300 


,671 





Table 5» Correlation Coefficients for 44,000 > 
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Figure 10. Graphs of simple linear equations applicable to T 3<i>^ oca 
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Figure 11. Nomogram for = - /JZ3 4- L 
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Figure 12. Nomogram for 








Mean 



Standard 

Deviation 



Zj 36^518 3,809 

^LMW 35,169 3,372 

^300 29,773 594 

T-fble 6o Means and Standard Deviations for 
44,000 > ac:>o' 

and ^ 



fhe regression equations and their standard errors are as 
£cl lows 



^ iZ& vS. t , 

^ S/ & ^ 

csi/ T^- ^ 



^ /^TT- 

>se. - ^/?o 



(iZ) 

(f3) 

( 1 4 -) 



The graphs of equations (12) and (13) are drawn in Figure 13, 

A noraogr«^ 1 D for equation (14) is shown in Figure 14« The value 
of the independent variable of equation (14) is limited by the 
fact that their coefficients are small, thus giving an almost 
horizontal planCo Equations (4) and (5) show that the values 
of b and c approach zero as the three linear correlation co* 
efficients become equal « Table (5) shov/s the near equality of 
the values of r resulting from this particular selection of 

d 3i 1 5' o 
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Figure 13. Graphs of simple linear equations applicable to 

2-j.^ So^i^oa'' 
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Figure 14. Nomogram for ^ -h , Jlo? 
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Conclusions ■and Recommendations 



The assuTption that Zj = is considered acceptable 

onJv where Z-j. lies between 36,000 and 40,000 feet. Above and 
below thi-3 lange, the error associated with this assumption 
increases rapidly. Within the continental' United States, it is 
not ontomnon for the LMW to vary from 20,000 feet below the 
trcpcpause in the southern states to 12,000 feet above the 
tropopause in the north. 

It is concluded, then, that the muLtiple regression equa- 
Cions (7) and (ll)j, which together are applicable for all 
values of will give the best first approximation to 
There v/ill undoubtedly be weather systems which will be 
associated with large deviations of from the plane of 

regressioHa The conditions under which these deviations will 
occur , howevecg has not been investigated in this paper. 

The regression plane associated with equation (14) is 
very nearly horizontal and as a result, the two independent 
variables have small coefficients and thus very little effect 
on the value of the dependent variable, Z^j^« Further investi- 
gation with other independent variables In the range of 
applicability of this equation was prevented by insufficient 
t Ime 0 

It was originally intended to evaluate the regression 
equations using independent data from the month of December 
1959, Hov;ever, because of the availability of an electronic 
computer at the United States Naval Postgraduate School, it 
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was felt that time would be utilized more advantageously by con 
tinuing to search for the most significant correlations among 
the variables under consideration. In addition, it is believed 
that the variety of weather systems occurring during the months 
of November and January, make these regression equations suf- 
ficiently general so that they may be used with confidence for 
independent samples of data. 

Were this investigation to continue, consideration of 
independent variables other than Z'l and 2390 would certainly be 
warrantedo Since the winds are dynamically associated with the 
thermal structure of the atmosphere, the temperature at the 
300-mb surface may well show better association with 2jj^^ than 
the height of the 300-mb surface. Some work has been done in 
attempting correlations with 1000-500 mb thickness values 
It would be logical to expect that thickness values for layers 
above 500 mb would give results as good or better. 

The further consideration of 2y^y, 2»j> and or any 

other associated variables should include serial correlations 
and graphical solutions of their non-linear relationships* 
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